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ABSTRACT 



The met lod of turbine performance prediction developed by 
Vavra and Eckert has been refined in this analysis to realize more 
of the potential of the three-dimensional calculating methods. 

Mach number and rotor tip clearance effects on blade outlet angles 
and loss coefficients have been localized rather than averaged over 
the blade height. An approximation for streamline curvature has 
been used. 

Performance curves were determined for two single stage axial- 
flow turbines located at the Propulsion Laboratory of the Naval 
Postgraduate School. Test results were available for one of the 
turbines. Agreement between predicted and experimental performance 
values was generally within 3 per cent. 



Thesis by Robert G. Harrison entitled: M An Analysis of Single-Stage Axial- 

Flow Turbine Performance Using Three-Dimensional Calculating Methods” 
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1. Introduction 

Turbines form an important part of propulsion systems. To opti- 
mize a design it is necessary to know the performance at off-design 
conditions as well as the performance at the design point. Since the 
testing of a prototype is very costly and time consuming, it is of 
great advantage to be able to predict the performance characteristics 
by means of theoretical methods. The more advanced a system is the 
more important it becomes to improve the accuracy of these methods. 
Hence it is necessary to base these methods on the fundamental laws of 
fluid dynamics rather than on ru .e-of-thumb approximations. The latter 
basis is possible only when data on previous designs are available. 

For new and advanced configurations, it will be necessary to apply re- 
fined methods which enable the designer and systems engineer to pre- 
dict the effect of proposed design changes. 

The principal equations that describe the flow properties in 
turbomachines are well known. These same equations are used as a basis 
for all proper "three-dimensional” calculating methods that have been 
cited in the technical literature. The methods differ however in the 
manner in which these equations are manipulated and applied. 

This thesis is concerned with the refinement of the three- 
dimensional method of analysis developed by Vavra and Eckert. Based on 
the physical dimensions of particular test turbines that are available 
at the Turbo-Propulsion Laboratory of the Naval Postgraduate School, 
performance curves were determined for these machines. Concurrently 
with the performance analysis, experimental tests were conducted on 
one of these turbines so that actual experimental results could be 
used to julge the accuracy of the proposed theoretical performance 
evaluation. 

In addition to his publications and classroom lectures which con- 
stituted the foundation for this analysis, Professor Vavra was very 
generous in providing guidance and counsel during the period of this 
work. For this I am greatly appreciative. I would also like to thank 
Lieutenants P. M. Commons and J. A. Messegee for making their experi- 
mental results available. 
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2. Basis for the Analysis 

The two conservation equations that must be satisfied to obtain 
a solution for the flow in a turbine are the equations of motion and 
continuity. These equations are satisfied at stations between blade 
rows. The method used is that given by Vavra. Vavra developed the 
equations of motion and continuity for absolute flows. The equa- 
tions in this form are readily useable for the position after the 

stator. Eckert later developed these equations for relative flows to 

2 

be used for rotor calculations. Eckert's conversion allows relative 
flow quantities to be handled directly without conversion to an 
absolute system. Eckert's approach also avoids iteration procedures 
to determine the total enthalpy after the rotor. 

The assumptions made for the developement of the equations 
used in the performance analysis are: 

1. An infinite number of blades in each row so that downstream 
effects are not felt upstream. 

2. Axisymmetric flow at the stations where the equations of 
motion are solved. 

3. Adiabatic and steady flow so that the total enthalpy along 
any given streamline is constant through the stator and the relative 
total enthalpy is constant through the rotor. 

4. All entropy changes are assumed to occur in the blade 
channels that are located ahead of the stations where the equations 
of motion are satisfied. Hence at the calculating stations the 
flow is assumed to be isentropic along particular streamlines. 

With the above assumptions, the equations of motion for 
absolute and relative flows, respecti /ely, are 

yH = v x (v* V WTvs (1) 

vH k = Wx(vxW + 2o5') +T VS (2) 



^Vavra, M. H., Aero-Thermodynamics and Flow In Turbomachines . 
New York, London: John Wiley and Sons, Inc. j 1960, Chapter 16. 

2 

Eckert, R. H., Performance Analysis and Initial Tests of a 
Transonic Turbine Test Rig (USNPGS Thesis, May 1966), pp. 149-155. 
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Relative total enthalpy can be written as 

U -U + wf U, £ _ h \N? L )| . l 

Hr" n, + -^ - a ~ n 

Equivalent e.rthalov i* defined as 
With ,1 _ Oi 



2-S 



w4 yl 
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H e = h, s + 



- Nr 



Similar to H f the equivalent enthalpy h £ is constant along a 



a " K z 

en 

streamline for the adopted assumptions 

j* 

•e. ''a ' a 



also be written as 

w/_ h + Wil 
+ _ _ n, + 2 _ 



Equivalent enthalpy can 



H e 



( 4 ) 



(5) 



The introduction of the equivalent enthalpy allows this quantity to 
be used for the rotor in a manner analogous to the way total enthalpy 
is used for the stator. 

For equations used in this analysis, the subscripts refer to: 

0 - station ahead of stator 

1 - station ahead of rotor 

2 - station after rotor 

is - isentroplc expansion from 

s - isentropic expansion from P ° 
th - the theoretical value 

Figure 1 is a temperature-entropy diagram showing the thermo- 
dynamic process along a particular streamline for a single stage 
turbine. In general, the fluid properties will vary from stream- 
line to streamline. The method in which the loss coefficients are 
applied is also indicated in Fig. 1, The loss coefficients are 
defined as 



-v, 1 

X For the Stator 

•th 

(*7 _ Wgfh ~ Wg For the Rotor 

' 7r ‘ w4h 

The coordinate system that will be used in the analysis is 
indicated in Fig. 2. This figure also shows the general layout of 
the tvpe turbine to which the prediction performance analysis is 
applied , 

igr: convention for the various angles that are needed in the 
analysis l indicated in Fig. 3. 




(6) 

(7) 



19 



The modification of Eq. 2 into a form that can be used for the 
analysis is given in Appendix A, Section 1. The appropriate form of 
Eq. 1 can be obtained from the modification of Eq. 2 if the angular 
velocity at is set equal to zero. Other differences of the final 
equations derived from Eqs. 1 and 2 are listed in Appendix A, 

Section 1. Equation 1 can then be written 



' d x, 



-i ijr 

, 11^7 'Jdx. 



Equation 2 becomes 

itMJl - l/W Sr\ zimf 4^.4. 

dXj. COS L ;J ^ t_ £ JdXiJ z ' dx a X £ 

4iUcss^. & 



SIM 






w^.Y, 

where: 



Wa 



w/ Y.* dx 






( 8 ) 



(9) 



V= & or ^ 

1 V W,„ 


for Eqs. 8 and 9 respectively 
(subscript m refers to mean streamline) 


K=S-0 


streamline curvature factor 


%r 


streamline displacements shown in Fig. 2 


aR 




, _ E, La 
L “ a 


(see Fig. 2) 



s * = -§_ 
b c r 

c, = a^j J 

The equation of continuity is used in its non-dimensional form 
by introducing a flow function ^ . Developement of ^ is given 

in Appendix A, Section 2. In differential form the equations for 
absolute and relative flows can then be expressed by 



and 



JM 



{si 


i(t: 


f- 


(-t) 


! *J 








■&)*] 



The differential element of area dA is 

J A= i c acW 



( 10 ) 



(ID 



( 12 ) 
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where : 



Z ° number of blades 
a ■ blade exit opening (see Fig. 46) 

£ - arc; restriction coefficient 

Since £ * s valid for isentropic flow only, the restriction 

factor 5 must be introduced to correct the actual flow area to an 

effective area which accounts for the restrictions due to the 

boundary layers on both sides of the flow channel. 

3 

The factor £' can be expressed by 
^21 - 1 

j = t - U H ***_ | + y_ 



,03) 



In Eq, 13, £ * is the boundary layer displacement thickness and 

H*** is the so-called energy parameter defined as the energy thick- 

4 . 

ness divided by the displacement thickness. The term repre- 

sents the loss that is assumed to occur from the inlet to the 
throat of the blade channel, where is the loss coefficient repre- 

senting all the losses across the row of blades. The profile loss 

4 

coefficient was used by Eckert to represent the loss prior to the 
blade throat. Percentage of the total loss due to profile losses 
will vary considerably depending on blade geometry, radial position, 
and the incidence of the flow on the leading edge of the blade. 

Since secondary flow and tip clearance effects result in losses in 
the blade channel, half the total loss coefficient provides a better 
average representation of the losses in the blade channel prior to 
the throat. Tl e basis for the developement of r and H lV ** as used 
is given in Appendix A, Sections 3 and 4. 

by multiplying and dividing by l /v , Eq. 12 is 



-< .-*v 



(U) 



"vavra, M, ii, , Problems of Fluid Mechanics in Radial Turbo- 
machines (Khode-Saint -Cenese, Belgium: Von Karmin Institute for 
Fluid Lynamics, 1965) VKI Course Note 55b, pp. G46-50. 

^Eckert, oj>, cit . , p. 44. 
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After integration Eqs. 10 and 11 become, respectively, 



H d * 



( 15 ) 



and 



* L . Xh 



( 16 ) 



The flowrate w can be computed from the conditions ahead of the 
stator. Then a reference flowrate is defined by 

w^Tto R I 

PJ5 j 



'AVcf = 



(17) 



where W is in square inches. Continuity will be satisfied for 
the stator by 

At 1 

| idX =W, 

** ■* STATOR 

Similarly for the rotor 

» A r 

cl , _ / I -&£ 



ref 



p t e 



■p- -a- |}dX 



= W 



ref 



(18) 



(19) 



- ROTOR 



The in luence of the leakage flow through the radial tip clearance 
has not been accounted for in Eq. 19. 

The element of area between the blade tips and the shroud is 

J P\ ~ atr rdr 

The flow through the tip clearance area is 

,*s 

= 2 r, 



(70) 



id i-c R 

a, 1 3 



■/ 



x 



I i JX 

the value 

and X 2 for the tip will be used in Eq. 21. 
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,p 

*" _EftF^NlCE 

Since the tip clearance is relatively small, the values of 

. T. 



i. 5. 

f ith thesi 

,(\lf sj «*), 



With these assumptions Eq. 19 can be expressed by 

_ 



di 

r TlP 2 a » 



R(To t 



- tCCL r ,,\ 

L a m2: VmlFOToR ' 

The assumptions used to arrive at Eq. 22 are obviously incorrect in 
two respects. First, the flow represented by | r is not 
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perpendicular to the tip clearance area. Second, the effective area 
represented by ^ is larger than that which probably occurs 
because of the relatively large boundary layers that exist on the 
shroud and blade tips. The exact behavior of the flow in the small 
region between the shroud and the blade tips is impossible to predict 
without further tests. However, it is felt that the flowrate through 
this space as represented in Eq. 22 is too large for the reasons just 
mentioned. Therefore a more accurate approximation of this flowrate 
will be obtained if the last term on the left side of Eq. 22 is 
divided by 2, yielding 

f -Pit., Jjo £ S cl X •+• TTYV 
v 1 




Wr 



vre-f 
r L «^Hr KA ] p 



J ROTOR 

The tip clearance flow included in Eq. 23 can be obtained by inte- 
gration, 



(23) 



TT fy 



; Atli I i S } It MJ 






k _ rrJefr 



{ Pel 'I Ttl § — )ti i 



(24) 



Introducing this expression into Eq. 23 gi ves 

f / *AL,lST a. r I j)( + : / 3aJj£ e l ) 1 

Uk «*. Ute «~S* dA l Pf. 1 T Ie At,pJ 



ROTOR 



- 

= K2^] E . roR (25) 

3. Technique for Obtaining Solution 

With the equations of motion and continuity in the forms given 
by Eqs . 8, 9, 18 and 25, a method has been developed to analyze 
single stage axial turbines, in particular, those available for test 
in the Turbine Test Rig of the Turbo-Propulsion Laboratory of the 
Naval Postgraduate School. The method of analysis predicts turbine 
performance for specified values of inlet total pressure, inlet 
total temperature, rotor speed and the ratio of total inlet to static 
discharge pressure /P^. This method is similar to that described 
by Eckert.'* However, °Eckert ’ s analysis neglected some effects which 



5 Ibid, Section 3. 
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have been accounted for In this developement. Some significant 
changes made in this method are listed below: 

1. Stator and rotor outlet angles for a particular radial 
location are computed using a calculated Mach number for that loca- 
tion rather than an assumed Mach number or the Mach number of the 
mean streamline. 

2. Variation of loss coefficients due to changes in blade geo- 
metry in the radial direction is accounted for. 

3. The Influence of rotor tip clearance on the rotor outlet 
angles and loss coefficients is concentrated near the tip of the 
blade rather than averaging these effects over the full blade height. 

4. Fourth order polynomials are used to better approximate 
the curves respresenting blade characteristics as a function of 
radius and the curves of rotor loss coefficients as a function of 
incidence , 

5. Streamline curvature effects have been accounted for in the 
solution of the equations of motion. 

In addition to the assumptions that were mentioned in Section 
2 for the developement of the particular form of the equation of 
motion, the conditions ahead of the stator are assumed to be uniform; 
that la, the total temperature, velocity, and entropy are assumed 
constant and the flow axial in direction. It is realized that com- 
pletely uniform conditions are difficult to obtain, but any other 
assumpt on would be extremely difficult to develop mathematically. 

Direct solution of the equations of motion is not possible since 
they are nonlinear in the dependent variable Y. Likewise, no direct 
method is possible to eatlsfy continuity. Solutione of these equa- 
tions must therefore be gained by making initial assumptions for the 
values of the axial velocities which muat ba improved by successive 
iterations until the aquations are aatisfiad. To account for stream- 
line curvature and slope, a complete solution of the flow through 
the stator and rotor must first be made by neglecting the effects 
of curvature In order to determine atraamlina location* , Then the 
iteration to account for these affects may prograaa. These require- 
ments make the use of a high-speed computer « necessity. 
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This analysis has been programmed for the IBM 360 computer 
using FORTRAN IV . The program is described in Appendix C. The 
following paragraphs set forth the procedural steps of the program. 
The equations are listed in general form without referring to 
specific streamline locations. In the interest of clarity, however, 
some relationships will be written in forms similar to those used 
in the program. For example, P } (Z)/P t ^[^(2)/!^] *->* wil1 repre- 

sent the isentropic relationship for the number 2 streamline. 

Five streamlines are utilized for the analysis with the 
number 1 streamline located at the hub and the number 5 streamline 
located at the tip as shown in Fig. 2. The number 3 streamline will 
be used as the mean streamline, and the radius of this streamline 
will be referred to as the mean radius. The radial locations of the 
streamlines ahead of the stator will be such that the mass flowrate 
between adjacent streamlines is 25 per cent of the total flowrate. 
Positions for the streamlines after the stator and after the rotor 
are initially assumed. The locations of streamlines 2, 3, and 4 
then vary during the solution as necessary so that the percentage of 
the total flowrate between adjacent streamlines does not change. 

This continuity requirement will be called streamline continuity. 

Besides the radii, sufficient input information must be used 
to effectively reflect the physical characteristics of the stator 
and rotor blading. Some of the physical properties are introduced 
directly; such as, the number of stator blades, the number of rotor 
blades, and the rotor tip clearance. The other quantities used which 
reflect blade characteristics are throat opening dimensions for the 
blade channels, discharge angles, rotor blade inlet angles, loss 
coefficients, and stalling incidences for the rotor. 

Throat opening dimension "a" is a function of radius. The 
best method for introducing this characteristic into the anaylsis 
is to enter the measured values of "a" together with the correspond- 
ing radii. Then, utilizing the method of least squares, a fourth 
order polynominal curve is fitted through these points. From the 
resulting polynomial, the value of "a" for any radius required by 
streamline continuity can be determined. 
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Discharge angles are predicted by using a combination of the 
methods of Vavra and Ainley,^ Outlet angles are first calculated 
using the formula which Vavra established from the experimental data 
of Beei.^ These angles are then corrected for tip clearance, blade 
curvature, and Mach number effects with the methods given by Ainley. 
Stator and rotor discharge angles are predicted at three radii; 
namely, the hub, mean radius, and tip. The method used for calcula- 
ting these angles is explained in Appendix B, Section 1, 

Values of stator gas outlet angles for the mean stream- 

line are determined ior Mach numbers M^ of 0.5, 0.7, 0,75, 0.8, and 
1.0. These values are represented by two parabolic curves of the 
form 

^ = «+ bM, J - c M " ( 26 ) 



The first curve is used for Mach numbers M^ from 0.5 to 0.75 and the 



for Mj between 0.75 and 



second established interim values of CX. , 

1,0. From these curves the flow angle oc, 
can be found for any Mach number M, . 

The flow angle oc. is also a lunction of radius r^ There- 
fore the changes of oc, cor the hub and tip with reference to the 

mean radius, called and respectively, must be used. 

The flow angle <X, 
a linear variation ot 



for the mean streamline 



ran then be determined for any r^ by assuming 



CK 



between the hub and the mean radius and 



between the mean and the ti^, With this assumption and using the 
Mach number in Eq. 76 corresponding to the number 2 streamline, 



the flow angle cx, 



CK j (:Z) = OC | 



for this streamline would be 

jt (?) -ji*‘ 



(27) 



^A in ley, D. 0, and Mothiason, G. C, R,, A Method of Performance 
Estimation for Axial-Flow Turbines . Aeronautical Research Council, 

In 6 M No. 2974, 1957. pp. 3-4. 

^Beer, R,, Aerodynamic Design and Estimated Performance of a 
Two-Stage Curt is Turbine for the Liquid Oxygen Turbopump of the M-l 
Engine. NASA CR 54764 (AGO 8800-12), Nov. 19, 1965. p. 29. 
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The superscript ** is used with r^ and c£ f in Eq. 27 to indicate 

the radius initially assumed for the mean streamline and the 
computed for that radius. Equation 27 can then be used throughout 
the analysis even though the radial location of the mean streamline 
may change due to streamline continuity requirements. A similar 
approach is used to establish the flow angle at the rotor dis- 

charge. 

The rotor blade inlet angles are measured from the manu- 

facturing drawings of the blade profiles. Using the values of 
for the hub, mean, and tip streamlines, a parabolic curve is deter- 
mined which gives as a function of X^. 

Loss coefficients and stalling incidences are predicted by 

8 

using the methods of Ainley. For the present method, the stalling 

incidence i is defined as that at which the loss coefficient is 
s 

twice the value of the minimum loss coefficient. Following 
Ainley' s methods, loss coefficients are computed as a function of 
the ratio of flow incidence to stalling incidence . Loss 

coefficients are also a function of blade geometry or radius. Since 
the stator has zero incidence, its loss coefficient is a function 
of the radius only. 

Stator and rotor loss coefficients and rotor stalling incidences 
are calculated for the hub, mean, and tip radial locations. Fotor 
loss coefficients £ are determined for values of ranging 

from -2.0 to 1.6. Curves of vs. are drawn for each of 

the three radial locations, and values of , along with the 

corresponding quantities , are used to determine fourth 

order polynomials which approximate these curves. A similar proce- 
dur is followed to obtain a fourth order polynomial representing 
stalling incidence i g as a function of radius . Sample calcula- 
tions for the prediction of stator and rotor loss coefficients and 
stalling incidences are contained in Appendix B, Section 2. 



8 



Ainley, oj>. cit . , pp. 



4-5. 
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Variation of loss coefficients with radial location is accounted 
for by assuming a linear variation of these quantities between the 
hub and the mean radius and between the mean and the tip. To demon- 
strate the procedure followed in computing rotor loss coefficients, 
the following example is given. For a particular incidence, the 
first step in the determination of for the number 2 streamline 

would be to calculate i g by using the radius r, of that streamline 
and the polynomial of the form i fl * i g (r^). Loss coefficients for 
the hub and mean radius would then be computed using the resulting 

in the polynomials, for these radii, of the form = %('h) * 
With the assumed linear variation ( 2 ) would be 



V2>= Vo 



r ( **- r,co. 




<? 8 ) 



The reason for using the superscript ** on the mean streamline 
values is the same at previously mentioned in connection with 
8q, 27, 

For the first approximation , the Mach number M q ahead of the 
stator is assumed,and the static properties and flowrate at station 
0 are found by 



T _ 

lo 


(?9) 


V 0 =f7iln7 


(30) 


-■ - 


(31) 


a- k 


(3?) 


A,- trtr T ‘-r/) 


(33) 


w 5 ft* 


(34) 


$A* = Wref 


(35) 



The reference flowrate W ^ will be used to check continuity at 
stations L and 2, 
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then 



The next step is to determine axial velocities after the 
stator that satisfy the equation of motion. Total enthalpy after 
the stator is constant by assumption, and streamline curvature 
effects are neglected at this stage of the analysis. With these 
conditions, Eq. 8 simplifies to 

This equation can be integrated to give 

In Y, z *C*Idx + /nC t ' 

where X q is arbitrary and ffjC is the constant of integration. 

Using the boundary condition Y « 1.0 at X ■ 1.0, can be 

found by f f 

O- fldt* lhC z on lnC z =-J^Jc/X 

*o x. 

in y, z -fidx- Jm=fm 

X 0 Xc / 

Eouation 38 must be expressed in a form that can be utilized 

in the computer. Expansion by infinite series yields ^ 

Y= e ifM = , + „ + "?+ Jg. + g+ - n 

The auantities contained in I, Eq. 36, must be evaluated 
before proceeding with the solution. For the first approximation, 
assumptions are made for the values of Y^, M^, and . The flow 
angles are then calculated by using Eqs. 26 and m 27 . After the 

value of d. has been determined for each streamline, L is 

or X | 

computed. This derivative and all others needed in the analysis are 
found by finite difference methods. 

Enthalpy is computed by 

/-/= Cpl t 0 

9 

and the entropy term is found by the method of Vavra. 



( 36 ) 



(37) 



(38) 



(39) 



(40) 



Vavra, M. H. , Aero-Thermodynamics and Flow in Turbomachines . 
New York, London: John Wiley and Sons, Inc,, 1960. pp. 445-447 . 
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I 



V , 2 C, 

c \ Hco s*ot| 

y* v-r«. 



(41) 



C,Hco^Xh%)_ 

The stator loss coefficients are determined by equations similar 
to Eq. 28. 

Each time a solution to Eq. 36 is found, the new values of 
are then used for the next iteration. After five iterations, o( , 

doC 

and ■ are recalculated using the new stator exit Mach numbers. 
The Mach number for each streamline is found by 




M = , — - 

n ’ 'i *■ 9 RT, 

With the new values for o( 



and 



doc, 

dX 



(42) 
(42a) 

(43) 

(44) 



five more iterations are 



made to determine the corrected values of Y^. 

The quantities represented by Eqs. 42-44 are recomputed after 
satisfying the equation of motion and additional quantities deter- 
mined by 

Vu,= \ Tfl/JOi.! 

T . - T - 
'H%- 1 1 0 |-^ s 

\V«-i 



p. s p 4 lr) 

The flowrate through the stator is computed next and compared 
with that required to satisfy continuity. Reference flowrate between 
the hub and each streamline is found from 

/ / 

50 tA = a ^ r ^J | WI dX 



(«) 

(46) 

(47) 



(48) 



where 









f 



(49) 
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The "a u and "a^" in Eq. 48 are found from the polynomial which 
represents throat opening as a function of radius. 

Before Eq, 48 is solved, the pressure ratio is compared with 
the critical pressure ratio. If the critical pressure ratio has 
been exceeded, the flow is choked at that radial location and the 

T w> 

critical pressure ratio is used in calculating 4> and ~ 

The fraction of the total flowrate passing between the hub and 
each streamline is computed by 

WIJX 

/ WI Jx 

Total flowrate, as found by the denominator of Eq, 50 multiplied 
by , is then compared with the reference flowrate. 

Overall continuity is satisfied if the difference is less than 

0 . 0002 . 



If required and computed flowrates are not within this toler- 
ance, the axial velocity for the mean streamline is adjusted by 

W--o^:p'jted 



v,. 



(mew) "Vawwl *- r; ) o.ooos? 

W is the required reference flowrate divided by and 

% req 

W , is the denominator of Eq. 50. 

computed ^ 

Solutions to the equation of motion and continuity are 
successively found until overall continuity is satisfied. 

The fractions of the total flowrate determined for each 
streamline by Lq. 50 are then compared with the corresponding flow- 
rate fractions ahead tf the stator, Streamline continuity is satis- 
fied if agreement is within 0.002 for each streamline. If stream- 
line continuity has not been satisfied, the streamlines in error are 
adjusted by 



X, r ,, = X 0 ,o *[w, r . r w 

Equation 52 applies only to streamlines 2, 3, and 4 since = 0 at 
the hub and = 1.0 at the tip. 



(50) 



(51) 



(52) 
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If streamline positions have been adjusted, new streamline 
radii are found by 






NE.W OLD 



( 53 ) 



and new values of X for each streamline are obtained by 

X = ^ (54) 

With the new streamline radii the equation of motion and overall 
continuity must be satisfied again. 

The following values are determined after streamline continuity 
is satisfied: 

rrRFNi 
U| (30)(l£) 

U a - !M 

' I 



W.^Vu.-U, 

,, -1/ Wo i \ 

a = (_j 

IA/ _ Va, 

T = T * + 

H £ =T tL c P 



Pt 



*6 ’ T, 



T tt 



V*-i 



(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 



Rotor blade inlet angles are determined for each streamline 
by means of the parabola which establishes /@ 0 as a function of X^« 
Incidence is found by 

i=/2, - /6 a (63) 
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Stalling incidences and rotor loss coefficients are then deter- 
mined using the procdure described in connection with Eq. 28. If 

the incidence ratio -d — is less than -2.0 or greater than 1.6, 

: is 

the value for 

(s 



the rotor is 



cl s* 

.a y 



ds 



is set equal to -2.0 or 1.6, respectively. 

One of the quantities necessary for the equation of motion after 

i *■ 

d Ss - This quantity is separated into two parts, 



J X - 



ds~j 

c\t. 



( 64 ) 



represents the entropy gradient due to changes of entropy 

I ' * 

across the stator and referred to station 2. The term - c 

d Aa 

represents the entropy gradient due to the different entropy changes 
through the rotor. The entropy increase through the rotor is com- 
puted by using the corresponding values of the rotor in Eq. 41; for 
examplejHg would be used instead of H. 

Neglecting streamline curvature and slope, Eq. 9 can be re- 
written 

A(LV) _ 

d Xo 






z LU U 



cl _ 4 LW S I N @ 2 - 

d X& " 'i z ' 

a Cl dHe /. C, He COS~-^ a\ ds? 

v" y / dx t 



Y/ Wa.^ y, 2 w*L. c |y, V ^ /dXa < 65 > 

The discharge angle and its derivative are found in the 

d Xa 

same manner as previously described for the stator discharge angles. 

The method for solving Eq. 65 is similar to that used for Eq. 36 
with the exception that iterations are carried out until correspond- 
ing values of Y change by less than 0.005, or unv.il 13 Iterations have 
been completed. The extra iterations are necessary because there 
may be a slower convergence of the values of Y at station 2. After 
satisfying the equation of motion, values at station 2 corresponding 
to those represented by Eqs. 42-47, are computed. It should be 
noted that where absolute velocity terms are used in Eqs. 42-47, the 
corresponding equations for the rotor will employ relative velocities. 
Therefore is the Mach number of the flow relative to the rotating 
rotor blade. 
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Overall continuity is checked at the rotor discharge using 
the same procedure utilized for the stator; however, the reference 
flowrate is increased according to Eq. 25 to account for the flow- 
rate between the blade tips and the surrounding shroud. 

After overall continuity has been satisfied at station 2, 
streamline continuity is checked. In addition to the functions 
performed for the check after the stator, there are certain quanti- 
ties that must be adjusted for streamline relocation. These are 




hu wrv/ = Hr, 



.IX 



r(x 



\ Airy./ j 



F OUD d X; 



Mtvs/ Xou v) 

* (X„r 






•X.LO) 



A solution exists when the equations of motionjand overall 
and streamline continuity have been satisfied; however, this solu- 
tion has neglected streamline curvature and streamline slope. To 
obtain a solution which accounts for streamline curvature effects, 
the terms in Eqs. 8 and 9 that contain or L R must be in- 

cluded when solving the equations of motion. 

The streamline displacements Sr and 0R, shown in Fig. 2» 
are computed for each streamline by 



Sr* r, - 



r 0 - r, 

2 



and 



ap * r 0 - f? 

Equation 70 is based on the assumption that the radii of any stream- 
lint; at stations 0 and 2 are not greatly different. Therefore the 
cosine of the angle between the line connecting these points end the 
axis of the machine is approximately equal to unity. The length L, 
also shown in Fig. 2, is taken to be half the distance from 0.1 inch 
ahead of the stator to 0.1 Inch after the rotor. 



( 66 ) 

(67) 

(68) 

(69) 

(70) 

(71) 
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After a solution is obtained which accounts for streamline 

Pro 



curvature effects, the resultant pressure ratio 



compared with the 



Pto 



?Z 



specified initially. Ideally 



is 

Pro 



for 



Pa P 2 

each streamline will be the same. However, computer solutions for 

this quantity may vary slightly from streamline to streamline. 

p “ 

Therefore a mass-flow-weighted value of — , called 

Pz. 

is found by 







(72) 



If the specified and f— differ by more than 0.0003, the 

Pa \ ?z)\N 

Mach number of the flow ahead of the stator is properly adjusted 
and another solution is found. 

P 

After the iterations for — have been completed, additional 
quantities are determined for each streamline , using 



AH = H,-H a - (U.Vo, - U^VoJ) 3 y 


(73) 


T - T - aH 

lt 0 c P 


(74) 


H 




W\(i) ‘ 


(75) 



Overall efficiency is then computed by 

L H _ Tt 0 ~ Tr a 

a hi* Tffc-Tk^ 

The ideal change in enthalpy A his is the isentropic enthalpy 

drop from the total inlet pressure F^ 0 to the static discharge 

pressure Equation 76 is used to compute the efficiency of a 



33 



single stage turbine because the kinetic energy leaving the rotor 

cannot be utilized. The efficiency defined by En. 76 will be 

referred to as total-static efficiency. 

Theoretical degree of reaction r* and head coefficient k. are 

is 

given by 



r 



* 



- Tzt% 



Tto ~ /zC5 



_ ZC 9 (Tj Q ~ 7zts)^H 
65 ~ = </,' 



( 77 ) 



(78) 



For turbine performance curves it is desirable to obtain the 

mass-flow-weighted values of efficiency ( )., head coefficient 

w 

(k. )., theoretical degree of reaction (r*)., horsepower (HP)., and 
is w w r W 

moment (>L) . . The last two nuantities are found from 



/jjp) _ (Atlhii JW 

~ •SSO 



and 



Mi - 



(h/P),;, (sso) 



CO 



The mass- flow-weighted A bl , called ( ^ H)^, as we ll as ( ^ )^» 

(kis)w’ anc * ( r *)^ 7 are com P uta d using enuations similar to Er?. 72. 
Referred values are obtained, following NASA practice. 
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For -1.4: 



Np^ ef - (* p )* 
Yef fe £ 


(81) 
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where: 



g = H 2 . - 7~ to - 

Tstd. 5 / 8.4 

( 85 ) 

C _ ^0 _ _^D 

d “ Pi-ra. A7.7 (86) 



4. MOD I and MOD II Turbines 

The method of analysis as presented was used to determine 
performance curves for the so-called MOD I and MOD II turbines. 

Both turbines are single stage axial-flow machines. Experimental 
tests were conducted on the MOD II turbine by Commons and Messegee 
and are described in Fefs. 4 and 6. The test results are plotted 
with appropriate predicted performance curves. 

The so-called MOD I turbine was designed for free-vortex flow 
and has highly twisted blades. Outer diameter of this turbine is 
9.898 inches. The hub diameters of the stator inlet and rotor dis- 
charge are 6.930 and 5.970 inches, respectively. The stator con- 
tains 13 blades and the rotor 22 blades. Blades of the MOD I tur- 
bine are generally thin. The stator and rotor profiles used to 
predict outlet angles and loss coefficients are shown in Fig. 4. 

The MOD II turbine is approximately the same size as the MOD I, 
but its blading is considerably different. The blades of the 
MOD II turbine are thick with blunt leading edges and constant 
profiles over the blade height. Outer diameters of the MOD II 
turbine stator and rotor are 9.701 and 9.836 inches, respectively. 
The stator has a hub diameter of 6.796 inches, and the hub diameter 
of the rotor is 6.598 inches. There are 19 stator blades and 18 
rotor blades. Stator and rotor blade profiles for the MOD II tur- 
bine are shown in Fig. 5. 

Throughout the remainder of this thesis the MOD I and MOD II 
turbines will be referred to simply as MOD I and MOD II. 

The minimum throat opening "a" of the blade channels is a very 
critical dimension. Slight variations of this quantity have a con- 
siderable effect on turbine performance. Since this quantity is so 
sensitive, values of "a" measured from the actual hardware were used 
for the analysis rather than those obtained from the manufacturing 
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drawings. Then errors due to manufacturing will not be a factor in 
comparison of predicted and experimental results. Figure 6 shows 
the throat openings "a" as a function of radius for the stator and 
rotor blading of both turbines. 

Predicted stator outlet angle as a function of radius r^ 

for the MOD I is plotted in Fig. 7. The assumed linear variation of 
between the hub and the mean radius and between the mean 
and the tip is readily apparent in this figure. Figure 8 shows the 
predicted variations of \ , with Mach number M^ for the MOD I. 
Figures 9 and 10 are the corresponding plots for the MOD II. 

Relative discharge angles / S 1 were computed for two radial tip 
clearances k for each turbine; namely, 0.020 and 0.033 inches for 
the MOD I, and 0.015 and 0.033 inches for the MOD II. The predicted 
flow angles are plotted in the same manner as previously 

described for the flow angles c\ . Figures 11 and 12 show 
as a function of radius r^ and as a function of Mach number M^, 
respectively, for the MOD I, where M^ refers to the Mach number 
of the flow relative to the rotor. Figures 13 and 14 show the 
corresponding plots for the MOD II. The predicted effect of radial 
tip clearance on the discharge angles can seen in Figs. 11 - 

14. 

Stator loss coefficients were computed for the radial 

locations corresponding to the hub, mean radius, and tip. Figure 
15 shows as a function of radius r^ for both the MOD I and 

MOD II. The straight lines in this figure between the values of 

at the hub and mean radius > and between the mean and the tip , 
reflect the assumed linear variation of in these regions. 

Variation of rotor blade inlet angle ^ \ with radius r^ is 
plotted for both turbines in Fig. 16. The difference between the 
untwisted and the free-vortex blades is easily noted in this plot. 
Also shown in this figure are curves representing the variation of 
stalling incidence i g with r^. The change of the MOD I rotor 
blade profiles with radius is reflected by a considerable variation 
of i g whereas just the opposite is true for the MOD II. 
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Curves for the MOD I showing predicted rotor loss coefficients 
as a function of incidence ratio -4 — for the hub, mean radius, 
and tip are shown in Fig. 17. Since the loss coefficients between the 
mean radius and the tip are dependent on tip clearance, there are two 
curves for the tip. One curve holds for the tip clearance of 0,020 
inches; the other is for the larger tip clearance of 0.033 inches. 

For negative incidence ratios between -0.5 and -2.0 the curves for 
the tip are estimations. This was necessary because the computations 
by the method shown in Appendix B, Section 2, gave unrealistically 
low values of as the flow inlet angle ^ approached -90°. 

The situation is more easily understood when it is noted that the 
blade inlet angle is -33.7° and the stalling incidence is 37° at 
the rotor blade tip. Figure 17 shows that the loss coefficients for 
the hub (r^«3.300 in.) are relatively large. The loss coefficients 
at the hub are larger for most incidence ratios than those at the 
tip for a tip clearance of 0.020 inches. The larger value of 
at the hub reflects the higher losses that are associated with an 
impulse type blade. It can be seen in Fig. 4 that the blade shape 
varies from a reaction type profile at the tip to an impulse type 
profile at the hub. 

Loss coefficients for the MOD II rotor are plotted in Fig. 18 
for tip clearances of 0.015 and 0.033 inches. The predicted 
similarity of the curves for the hub and mean radius is to be ex- 
pected since the blading differs only in solidity. Although the 
blade profile is the same at all radii, the losses due to tip 
clearance result in larger loss coefficients for the tip profile. 

For convenience, the blade properties used for calculating the 
MOD II rotor loss coefficients were those at the hub, mean and tip 
radii of the rotor discharge. Since the annulus area at the rotor 
discharge is larger than the annulus area at the stator discharge 
and since the flow incidence is a significant parameter for the 
rotor loss coefficients, it would have been more appropriate to use 
the blade characteristics at the hub, mean, and tip radii of the 
rotor inlet. However, the error is insignificant because the blade 
profile does not change along the radius. 
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It may be noted fiom Fig. 4 that the minimum radius, for which 
a blade profile is given, is 3.597 inches whereas the radius at 
the hub of the MOD I stator ex^ t is 3.300 inches. The outlet angle 
and loss coefficient for the hub were found by extrapolation, using 
the values computed for the radii of 4.125 and 3.597 inches and 
assuming a linear variation of these quantities with radius. The 
relative outlet angle for the hub of the MOD I rotor was found in 
a similar manner. 

Performance curves for the two turbines were determined for 
the rotor tip clearances mentioned earlier. The total inle t to 
static discharge pressure ratios investigated were 1.30, 1.40, 

1.50, and 1.60, with the exception that pressure ratios of 1.31 and 
1.51 were used for the MOD II with 0.033 inch rotor tip clearance. 
These pressure ratios agree more closely with those experimentally 
investigated by Commons and Messegee. 

The axial distances L used for the determination of the 
curvatures depend on the axial clearance between the stator and 
rotor as well as on the blade geometries. Axial clearances of 0.4 
and 1.0 inches were used for the analyses of the MOD I and MOD II, 
respectively. 

Curves representing the performance of the MOD I are plotted 
in Figs. 19 through 26. Performance values plotted are mass- flow- 
weighted values unless stated otherwise. Figure 19 shows referred 
flowrate as a function of referred RPM. The increase in flowrate 
due to an increase in rotor tip clearance can be seen in this 
figure. Although the flowrate is greater for the larger tip 
clearance , the toi que developed is greater at the smaller tip 
clearance. The decrease in torque for the larger tip clearance 
results from the increased losses and decrease in turning angle of 
the flow through the rotor near the tip. The predicted effect of 
tip clearance on torque is shown in Fig. 20 where referred moment 
is plotted versus referred RPM. 

The variation of total-static efficiency with referred RPM for 
the two tip clearances can be seen in Figs. 21 and 22. The referred 
RPM at which maximum efficiency occurs increases when the total 
inlet to static discharge pressure ratio is increased. Blade losses 



40 



and the kinetic energy of the flow leaving the rotor affect the total- 
static efficiency. As pressure ratio is increased the absolute veloc- 
ity leaving the stator and the relative velocity leaving the rotor 
increase. Therefore the peripheral speed of the rotor must increase 
to obtain conditions where the absolute velocity leaving the rotor 
is in an axial direction and where the relative flow ahead of the 
rotor has zero incidence. The RFM where the flow has zero incidence 
on the rotor will not necessarily be that at which the absolute 
velocity leaving the rotor is axial. At any RPM, flow incidence and 
absolute discharge angles vary from streamline to streamline, and the 
above statements refer to mass- flow-weighted values. 

It may be noted in Figs. 21 and 22 that the peak total-static 
efficiency decreases somewhat as the pressure ratio increases. At 
higher pressure ratios the ratio of kinetic energy leaving the rotor 
to the work done on the rotor increases. Since the kinetic energy 
leaving the rotor is lost energy for a single stage turbine, the 
total-static efficiency declines. The effects of pressure ratio and 
tip clearance on efficiency can be seen in Fig. 23 where total-static 
efficiency is plotted as a function of the isentropic head coefficient. 

The variation of referred power with referred RFM can be seen in 
Fig. 24. Peak power does not occur at the same referred RPM at 
which peak efficiency occurs. The peak referred power occurs at the 
referred RPM where the product of total-static efficiency and referred 
flowrate is greatest. 

Theoretical degree of reaction is plotted as a function of isen- 
tropic head coefficient in Fig. 25. It may be noted that the theo- 
retical degree of reaction increases with increasing pressure ratio and 
decreases with increasing tip clearance for any given isentropic head 
coefficient. The predicted effects are considerably different from 
the results of the radial turbine tests conducted by Riley. ^ 

Riley found that theoretical degree of reaction was independent of 
pressure ratio and axial clearance for radial turbines. 



Riley, M. W., The Effect of Axial Clearance on the 
Performance of a Dual Discharge Radial Turbine (USNPG Thesis, 
L/ecember 1966), p. 70. 
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Performance values for each streamline are obtained from the 
computer solution. However, plots utilizing values for each stream- 
line would be difficult to analyze. The deviation of the hub and 
tip values from that of the mass-flow-weighted average may be seen 
in Flo. 26. In this figure, hub, tip, and mass -flow-weigh ted values 
of theoretical degree of reaction are plotted as functions of 
referred RPM for a total inlet to static discharge pressure ratio of 
1.40. 

Performance curves for the MOD II corresponding to those pre- 
sented for the MOD I are plotted in Figs. 27 through 39. Additional 
plots have been used for the MOD II because of the inclusion of 
experimental results. An axial clearance of 1.0 inches was used for 
the theoretical prediction. Therefore, only experimental data for 
that axial clearance are shown. 

Comments made concerning the performance curves of the MOD I 
are applicable to the MOD II performance curves also. Differences 
in the performance of the two turbines will be discussed later. 

Plots of the variation of referred flowrate with referred RPM 
are shown in Figs. 27 and 28 for two rotor tip clearances. The 
quantitative values as well as the curve shapes agree well with the 
experimental data. The maximum difference between predicted and 
experimental referred flowrates occurs at a pressure ratio of 1.51 
for a tip clearance of 0.033 inches. There are two experimental 
points for this pressure ratio and tip clearance that differ from 
the predicted curve by about 2 per cent. The predicted and experi- 
mental values for all pressure ratios and tip clearances have an 
average difference of less than 1 per cent. 

Figures 29 and 30 show curves of referred moment versus 
referred RPM. The trends expressed by the predicted curves ere in 
excellent agreement with the experimental data. Although the quanti 
tative agreement between theoretical and experimental values is very 
good for three of the curves, the experimental torque is generally 
lower than the predicted torque. The average difference between pre 
dieted and experimental values is about 3 per cent. 



Figures 31 through 34 show total-static efficiency as a func- 
tion of referred RPM. The shapes of the predicted curves generally 
agree well with the experimental results. However, there is an indica- 
tion that experimental efficiencies decrease more rapidly at high RPM 
than is predicted by the theoretical curves. In the high RPM region 
the upper part of the rotor blade has a large negative flow incidence. 
In the prediction analysis when the incidence ratio had a value 
less than -2.0, the value of -2.0 was used for computing loss coeffi- 
cients. This limitation may be the reason that the predicted effi- 
ciencies in the high RPM range do not decrease as rapidly as the test 
data indicate. 

The quantitative agreement between predicted and experimental 
efficiencies varies considerably between different pressure ratios 
and tip clearances. There are two data points at a pressure ratio 
of 1.31 and tip clearance of 0.033 inches where the experimental 
efficiencies are over five points below the predicted values. At a 
pressure ratio of 1.50 and a tip clearance of 0.015 inches the 
average difference between predicted and experimental efficiencies 
is 1.5 points. Giving equal weight to all experimental values the 
average difference between experimental and predicted values is 
2.6 points. 

The calculations gave a decrease in efficiency by about two 
points as the tip clearance was increased from 0.015 to 0.033 
inches. The decrease in experimental efficiencies for the increased 
tip clearance varied with the different pressure ratios. However, 
the average decrease in efficiency is close to the predicted decrease. 

Total-static efficiency as a function of isentropic head co- 
efficient k^ g is shown in Fig. 35 for pressure ratios of 1.40 and 
1.60. Figure 38 shows degree of reaction r* versus k^ g . No experi- 
mental data are plotted in these figures because experimental values 
of mass- flow-weighted r* and were not available. 

Plots of referred power as a function of referred RPM are shown 
in Figs. 36 and 37. The comments made earlier concerning the referred 
moment plots apply to these curves also^ since the turbine power is 
proportional to the product of torque and RPM. 
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Hub, tip, and mass- flow-weighted values of theoretical degree 
of reaction are plotted in Fig. 39 for a pressure ratio of 1.40. 

The experimental values for the hub and tip cr 2 also plotted in this 
figure. The trend of the experimental data is in close agreement 
with the predicted trend; however, the quantitative agreement is poor, 
especially at the hub. The decrease in degree of reaction, as the tip 
clearance is increased ; is considerably greater than predicted. Both 
predicted and experimental results indicate that an axial flow tur- 
bine diffeis from a radial turbine inasmuch as the theoretical degree 
of reaction changes if rotor tip clearance is changed. 

Figure 40 shows the absolute flow angles and CX - as a 

function of radius for the stator and rotor outlets, respectively, 
of the MOD II. Experimental data were obtained for the same opera- 
ting conditions; namely, P /P =1.40,k=0. 015 inches, RPM/^ =13,394. 

t 0 £ 

The experimental and predicted stator outlet angles have an average 
difference of less than half a degree over the blade height. Pre- 
dicted and experimental agreement for the absolute flow angles at the 
rotor outlet is not as close. If the experimental point near the hub 
is neglected, the average difference between predicted and experimen- 
tal values of a is about 3.5 degrees. The experimental flow 
angles c*. are larger than the predicted angles over most of the 
blade height, which explains why the predicted torque was generally 
larger than the measured torque. 

Velocity distributions at the stator and rotor exits of the 
MOD II are shown in Fig. 41 for the operating conditions listed in 
the preceding paragraph. The experimentally found stator discharge 
velocities varied only slightly more from hub to tip than predicted. 

The average difference between predicted and experimental stator exit 
velocities was about 1 1/2 per cent. 

The experimentally determined absolute rotor exit velocities, 
which are shown in Fig. 41, are considerably larger than the pre- 
dicted velocities except at the hub and tip. Although experimental 
points could be taken only at distances greater than 0.16 inches from 
the inner diameter of the shroud, the trend of the data points indicates 
that the experimental velocities are less than predicted at the tip. 
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Kear the hub the experimental velocities were less than predicted. 

The decrease in experimental velocities near the hub and tip is 
the result of separated flow in these regions. The reader is 
referred to Kef. 6 for additional information concerning the experi- 
mental data and for photographs of the rotor showing indications of 
separation. 

Experimental relative rotor discharge angles /£? were deter- 
mined for the known rotor speed from the measured values of absolute 
rotor outlet angles and velocities, which are plotted in Figs. 40 
and 41, respectively. The experimental and predicted values of /S. 
are shown in Fig. 42. The magnitudes of the experimental values 
of near the hub are larger than predicted. Over approximately 

the outer three fourths of the rotor blade height the magnitudes of 
the experimental values of are less than predicted. At a 

radius 6f 4.4 inches, for instance, the experimental value of 
is about -57 degrees whereas the predicted value is about -67 
degrees. The average difference between experimental and predicted 
values of .r/ is about 6.5 degrees. The lower values of ^ pre- 
dicted near the hub,provide some compensation in the overall performance 
for the high values of predicted for the outer part of the blade. 

The MOL II was not designed to achieve the highest possible 
efficiency. The objectives of the design were to investigate the 
effects of b unt untwisted bladings of constant profile. This type 
of blading has the following advantages: 

1. Blade cooling passages are easily accommodated. 

2. A wider operating RRM range is posible for a specified tur- 
bine efficiency variation. 

3. Constant profile blades can be manufactured more economi- 
cally than twisted blades, especially if exotic high strength materials 
are needed for elevated temperatures. 

The following discussion concerning comparison of the MOD I and 
MOD II is based entirely on the predicted performance of these tur- 
bines using the methods described in this thesis. The reader should 
be aware that Ainley's methods for predicting loss coefficients were 
not developed for use with rotor blades having blunt leading edges. 
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Therefore the accuracy of the rotor loss coefficients predicted 
for the MOD II is probably not as good as for the MOD I. 

Figure 43 shows the predicted axial velocity ratios as a 
function of radius ratio at the stator exit and at the rotor exit 
for the MOD I and the MOD II. The curves are for the conditions 
that occur at the maximum efficiency of each turbine. The axial 
velocities of the MOD I are nearly constant over the blade height at 
the stator exit and at the rotor exit. The blading of this turbine 
was originally designed for free-vortex flow by assun ing uniform 
axial velocity components from hub to tip. From Fig. 43 it can be 
seen that this condition is only approximately satisfied since loss 
variations and curvature influences produce slight variations from 
the assumed distribution. Axial velocities for the MOD II decrease 
from the hub to the tip at the stator exit, and increase from the hub 
to the tip at the rotor exit, since the bladings of this turbine have 
constant profiles. The variation of axial velocities for the MOD II 

is more pronounced at the rotor exit than at the stator exit. At 

Was 

the rotor exit ~r varies from 0.735 at the hub to 1.195 at the 

tip. 

The maximum predicted total-static efficiency of the MOD I at 
a tip clearance of 0.033 inches is about 84 per cent. The MOD II 
has a maximum predicted efficiency at that tip clearance of about 80 
per cent. The difference in predicted peak efficiencies is due to 
different factors. The stator loss coefficients are highe for the 
MOD I than for the MOD II, whereas the loss coefficients for the 
MOD II rotor at zero incidence are somewhat larger than those for 
the MOD I at zero incidence. An accurate comparison of loss co- 
efficients is more involved for the rotor than for the stator. 
Because of its tw.sted blades ^the MOD I rotor has nearly zero inci- 
dence at all points along the blade height in the vicinity of the 
optimum efficiency. Only at one radius of the MOD II rotor blade, 
however, is the incidence zero. The incidence angles at larger 
radii are negative and at smaller radii they are positive. There- 
fore part of the blade is always operating at a loss coefficient 
larger than the minimum. Of probably even greater significance 
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for the efficiency decrease is the kinetic energy of the gas leaving 
the turbine. Minimum kinetic energy is lost when the flow is axial 
in direction. The design of the MOD I is such that at all points 
along the blade height the absolute velocity at the rotor exit is 
nearly axial in the vicinity of the optimum efficiency. The MOD II 
has only a small radial portion of the rotor blade where the dis- 
charge angle 0(^ is zero. At radii greater than the radius where 
OC 2 . is zero, the absolute flow discharge angle is positive and at 
smaller radii it is negative. Therefore the kinetic energy leaving 
the rotor of the MOD H is ,reater than that of the MOD I when both 
turbines are operating at peak efficiency. 

A comparison of the total-static efficiency versus isentropic 
head coefficient curves for the MOD I (Fig. 23) and MOD II (Fig. 35), 
shows a larger variation of efficiency with head coefficient for the 

p 

MOD I. For example, at k=0.033 inches and =1 . 60, the MOD I 

'Z 

efficiency decreases about 11.5 points as the head coefficient in- 
creases from 3 to 7 whereas the MOD II efficiency decreases only about 
7 points for the same change in head coefficient. This difference in 
efficiency variation results from the effects of blade twist mentioned 
earlier. 

To investigate the influence of the streamline curvatures on the 
predicted performance of the MOD I and MOD II the presented flow 
equations were solved for an axial length of the bladings L of 9 x 10”* 
inches. Increasing L to this value has the same effect as neglecting 
streamline curvature effects as can be seen from Eqs. 124 and 125. 
Differences between performance values neglecting streamline curvature, 
and those where curvature effects were accounted for^were less than 

0.2 per cent for both turbines. 

5. Conclusions and Recommendations 

Effects of streamline curvature were found to be insignificant 
for the MOD I and MOD II. The small effect on predicted turbine per- 
formance due to streamline curvature indicates that the exact value 
assumed for the curvature factor K will not greatly affect predicted 
performance. The results of this analysis also indicate that Vavra's 
method of approximating streamline curvature is of sufficient accuracy 
for methods of analysis where the flow equations are satisfied at 
stations between blade rows. 
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It is recommended that the MOD I be tested. Results of that 
investigation would provide additional information concerning the 
general applicability and accuracy of the method of analysis pro- 
posed in this thesis. 

Predicted and experimental flowrates for the MOD II were in 
close agreement. Since the restriction factor ^ is a significant 
factor in the predicted flowrates, the experimental results verify 
the validity of the theory used in the developement of J* 

Experimental results for the MOD II showed that the predicted 
rotor torque and turbine efficiency were generally 2 to 3 per cent 
too high. The angles and velocities predicted for the flow at the 
stator outlet were in excellent agreement with the experimental re- 
sults. Therefore the high values predicted for rotor torque and tur- 
bine efficiency must be tied to the rotor solution. The measured 
values of outlet angles and velocities for the rotor discharge also 
indicate that the predicted rotor solution is not completely correct. 

Exactly what parameters are in error in the predicted rotor 
solution, and to what extent, cannot be stated without more experi- 
mental data. More traverses should be taken at the stator and rotor 
outlets using recently calibrated probes. With this informat ion, loss 
coefficients and relative discharge angles for the rotor could be 
determined on a streamline basis. Based on the experimental results 
included in this thesis, it is suspected that the magnitudes of the 
predicted discharge angles 0 ^ are to ° large. If the predicted 
turning angles of the flow through the rotor were less, the predicted 
torque and efficiency would decrease. Also, the separated flow at 
the hub and tip at the discharge of the MOD II rotor indicates that 
predicted loss coefficients should be higher for the streamlines at 
these locations. 

To accurately predict turbine performance at off-design as well 
as design conditions, a streamline analysis is necessary. The experi- 
mental verification obtained for the proposed method indicates that 
the method possesses much potential and that additional developement 
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is warranted. For future developement work on this type of analysi 
it is recommended that the following changes be considered: 

1. Use seven streamlines instead of five. 

2. Calculate discharge angles * using the methods of this 
thesis but neglecting blade curvature effects. That would decrease 
the angles by 2 to 3 degrees and be in closer agreement with experi 
mental results. 

3. Introduce a multiplying factor for the rotor loss coeffici 
ents which would more accurately account for the separated flow at 
the hub and tip but would not greatly change the average rotor loss 
coefficient. For example, rotor loss coefficients would first be 
calculated using the methods of this thesis, and then the loss co- 
efficients for streamlines 1, 2, 6, and 7 would be multiplied by, 
say, 1.4 and those for streamlines 3, 4, and 5 would be multiplied 
by 0.6 or 0.7. 

4. Calculate rotor loss coefficients for an range of 

-3.0 to 2.0 instead of -2.0 to 1.6. The increased ~ range would 
improve performance prediction for turbines with untwisted rotor 
blades at off-design conditions. 
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TURBINE THERMODYNAMIC PROCESS 
FIG. I 
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VELOCITY DIAGRAMS 
FIG. 3 
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STATOR a ROTOR THROAT OPENINGS VS RADIAL POSITION (MODI a MOD II) 

FIG. 6 
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VARIATION OF STATOR OUTLET ANGLE WITH RADIUS (MODI) 

FIG. 7 
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VARIATION OF STATOR OUTLET ANGLE WITH MACH NO. (MODE) 

FIG. 10 
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VARIATION OF ROTOR OUTLET ANGLE WITH MACH NO. (MODI) 

FIG. 12 
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VARIATION OF ROTOR OUTLET ANGLE WITH RADIUS (MOD II ) 

FIG. 13 




VARIATION OF ROTOR OUTLET ANGLE WITH MACH NO. (MODI!) 

FIG. 14 
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STATOR LOSS COEFFICIENTS AS FUNCTION OF RADIUS (MOD I 8 II) 

FIG. 15 
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FIG. 16 




62 



ROTOR UOSS COEFFICIENTS AS FUNCTION OF INCIDENCE RATIO (MOD I ) 

FIG. 17 
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ROTOR LOSS COEFFICIENTS AS FUNCTION OF INCIDENCE RATIO (MOD II) 

FIG. 18 
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THEORETICAL DEGREE OF REACTION VS REFERRED RPM 
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APPENDIX A 



DEVELOPEMENT OF EQUATIONS 

1. Equation of Motion^ 

The equation of motion for relative flow is 

V H K = w X (v v Vv 1 2 Co) + T vs 

In cylindrical coordinates: 

~ ( v Tr 



x\ ; = 
r 



r 

- w„ 



* 7, )\ ^ 

^ . r i - i. 'b r 

~ ( v 

+ - r ' v ■ i r [ J p j d 5 

W x 2 uj = i r - dW;») - (2 w VY r ) 






[ 5VV j(rv\ ■':! / 

[ 5 © - 2 J j 






-3 



T , 



■ J_ 

Q r ao 



** X 1L + T ^-- N ) 

as ^ a* L r *r ) 



Equating Eq. 87 to the components in Eqs. 88-90, 
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1 
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r ae 
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( 2 ) 

(87) 



(88) 

(89) 

(90) 

(91) 

(92) 

(93) 
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Eckert, op. cit . , pp. 



149-155. 
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For assumed axisymmetric flow, ^ 

0= _y^>^_.vv r 



o. 



<' hi 

r 



^ r 



7. a Hr w aWr \ A/ SWh Wj ^rWJ) iT lS 
'^•7F' W ^T? ~ Wf FF r ~~ 1 



a* 



^2 



— ‘ W* ^-r + 2 uj\N 0 



■ • _ Wj ^cVWj ) _ \ a/ 

r * " r ^ " /VA 

Equation 94 can be written 

_ Wr 36-Wj) , f Wr. 

5? " ~ W ( ~ * W* 

Replacing — __1_ in Eq. 95 by its equivalent in Eq. 97, 

5 Wp \ a / 3'Af ... 3W* . Wj / Wr 3(rAj\\ _ . . Wj IVr . T 

~ = VVr tt ■ v " ~ * _ ft v«.>, * r r-^ w. 1 s* 

or 

“ ill. ,- - + Wr - Wr — -2 £ tT^- 

fVfe ;f r 32 v r 3f ' JJf. 3.- 

Multiplying Eq. 96 by and Eq. 99 by W^, 

. Mr „ vy-w- ^kx:l . w„ w r --Wr , w* w r — - -^ w, w r *w f T^ 

r o r V* ^ f 3 * « i 

±pL r _VvVli£ ^L/j^ 4. yV,. \Y r -WftWr -2-0 IA/WW *V/ fi T 

~ j?- r - ^ *r ^2 « f 



r S 

Adding Eqs. 100 and 101, 

Wr ; ' 



^2: 






For adiabatic flow is constant along a streamline, therefore, 

W dt ■7W P = c 



For axisymmetric flow, 

W J A -v H f - ~ = 



SWe. 3 Hr 

A + VVy- ~^TT 



From Eq. 104, 

ri r_ _ _ Wr ^ Hw 
- 2- Wa * r 

Substituting Eq. 105 into Eq. 102, 



Wr 7 ^7 7 t?) =t 7^ +w > 



^2 



VW 



(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

( 100 ) 

(101) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 
(107) 



90 



Using Eqs . 105 and 107, Eq. 99 becomes 

' A'r V - . v • f A i ^\ki ^Wr M ^ V W" i j> L . i Wj l > Wf T d S 

\f\p a r ~ rW/< r>r r ^ r Vy VV/i 

Multiplying Eq. 108 by f - — ^ gives 

' / . t J 



^ Hp _ V\y oin/^l _ yy _ 4 - vVn + 2 o Wj - 



T- 



<3 < 



>r 









(108) 



(109) 



Eq. 109 is identical with Eq. 96, therefore Eq. 109 is the equation 

f — | > Eq . 109 becomes 



that must be solved. With f VVp 



% (V v 

6f 



■K 



y'r' * + 






,-rlL - 0 

dr ’ 



( 110 ) 



Rewriting Eq. 110 for the rotor discharge and using the substitutions 



£ Cp 



The last term of Eq. Ill can be expressed by 

J_: -v..-vv (i )j=f p [ 2 H E -W^-VVpi y + ^)] 



With 



(-“ T ' ) and (' 'TAW h - 



Eq. 112 becomes 









Using Eq. 113, and after rearranging, Eq. Ill becomes 

*(A>. y ^ yA/ ^Sa VVW. Ayy 

c r 2 ~ y : F ^ <y 












Cp 






Noting that 

*(U?) v 



- =• ^ 2 



Eq. 114 is 

' i/.* . 



A v i ( - V V r. 2 + ^ j 

'' TF ^»y ~ f 

c> He ^ _ 

^ f * **" * 



Wd. T 4 ^ Wj- 



^ f 7 



v ‘ r - -i-l- -Wj* 1 — * - o 

C V- ' ~ 



;w.. 



(in) 



( 112 ) 



(113) 



(114) 



(115) 



(116) 
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/ \ 

This equation is made non-dimensional by multiplying by ( — )) 

where the subscript m refers to the mean streamline. In the sub- 
sequent derivations the subscript 2 will be omitted. Then 

(£_ W; \ 

Wj Ww. c ^ \Nt-.J 






M _ Mr 

2 M* £ 2 W^Y c p c O- r 



Ml ^ r 
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f YM 



Ml 



JTW 2 Hi 
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2He W* 1 

mi wi I 



3S 

Mi I a r 



= O 



Equation 117 then becomes 

C V ) _ , X ? V. ~ r A 2, oi" 

^ x / A/' cr ^ x 
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sx ° 



The fourth term in Eq. 121 can be expanded and rearranged as 
2 



yI M[lM Y --2 V^[xyYF -XTA^|1 + Yt^] 






= 2 Y T/M* /Y +2 Y T An /f 

> d A ' 
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2Y 



z „,ch_ 
-OX 



= T A H ‘ : 



v(Y 2 ) 



(117) 

(118) 

(119) 

( 120 ) 



( 121 ) 



(122) 
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Equation 122 then is 



^y 2) / v a y 2 r *s» y » ml Hi , 

0 + TflN, A) 2 ^ ccz 2 ^ 3X X T N ' 



, 4lUYrflM|g 2 ^vtc ^ iLU , 

kVfll ^x W/&. 



2 I-Ie 



.w A ; 



■Y Z rm z /S 
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Multiplying by ?-) > and with |+TAM^= ; ^ 7 ^ * 

\A^ ^2 "co^A ^>X /+* ^ 



J-*£f>W 
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4l U Sii^C os/g Zco%£ l 2 ^ Cc^ Z^ 

W*~Y wJy 4 " W£Y‘W VlCY" P 



— =0 



The terns and 



C 



■ 2 V^ 



5Wr Y 



( 122 a) 



w A 

streamline curvature and can be approximated by 



(123) 



represent the effects due to 



JL ^ i£ • Ks-s 

W A »* - * L* 1 



cos 2 X = 



. ‘Z /A^r\Z 

L + (^ 2 ") 

If the streamline curvature in meridional planes is zero. 



(124) 

(125) 



the terms 



represented by Eqs. 124 and 125 will take on values of zero and one, 
respectively. With $r positive as shown in Fig. 2 > the plus sign 
is used for K at station 2 and the minus sign at station 1. The 
streamline slope as represented by Eq. 125 is the same at stations 
1 and 2 since the streamline pattern is assumed to repeat itself after 
station 2 . 

Using Eqs. 124 and 125, Eq. 114 then becomes 






feKC. 



SKW- 



-Z T AM^f--f-S, N z /3 



W flw Y 



ZJXdJcos^ t a Co%& J|4e 2 lie eos^ 

vvtir 2 - w^r* jx 1 vo* ‘ Ti 



diT 

clX 



(126) 
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For the enthalpy terms to be In terms of BTU/lb m> they must be 
divided by 2gJ to keep the equation dimensionless. With 
C x - 2g J, Eq. 126 is 

4U^siM(6^/g zU^Ucof/8 C cos 2 /^ JHe fcUe*>»V 

W^iY* Wa^V 2 <IX [Wa"Y 2 f J d * 



4&,Y 2 )_ C6 a 



Noting the corresponding terms for the absolute flow of the 
stator, and with U=0, the equation for station 1 is 



dMl) 

dX 



co^a.If-Ks/'v 



Sr\ ( L 2 *-(i 

' L 2 / A l 2 - 
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- a TAMO*. 



c/s *■ 

4X 



- si mV + 



C, CoS ; oL 



sLH. 

dX 



C, HcosV 

Y 2 V^_ 



- S|N Z 0<- 



ds* 

dX 



2. Flow Function £ 

Flowrate in lb^/sec for the expansion process shown in 
Fig. 44 can be written 

w = 



(PLENUM) 




Expansion From Plenum 
Fig. 44 



(127)* 



(128) 



(129) 
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Assuming an isentropic expansion. 




then using the thermodynamic relations 

2 ijW P =T r T Cp= T TP 

the discharge velocity can be expressed by 




The density ^ is found with the isentropic assumption by 

JL . 

or 

i c \j - _ p 1 1 p t 

/’“M A l p t) KT r 

Using Eqs. 133 and 135, Eq. 129 becomes 

FL I PV 



= *k{-kY 



This can be rearranged to give 
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(130) 



(131) 



(132) 



(133) 



(134) 

(135) 



(136) 



(137) 



(138) 
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3. Restriction Factor 



12 



1 

The expansion process through a blade row is shown in Fig. 45 
(a) . The letter u represents the velocity at some point within the 
boundary layer and represents the velocity for an isentropic 
expansion. The various dimensions, angles, and velocities for the 
throat location are shown in Fig. 45(b). 








(fii) 



s — 
(>) 



Conditions at Exit of Blade Row 
Fig. 45 



The mass flowrate per unit blade height can be written 



where: 



- =A ^ s ' fe] + £ / V 

: ° 

* density at P^ for isentropic expansion (slugs/ft^) 

* density corresponding to the velocity u (slugs/ft^) 
<*d* - blade angle 

S ■ boundary layer thickness 



(139) 



^Vavra, M. H., Problems of Fluid Mechanics in Radial Tutbo- 
machines (Rhode-Saint-Genese, Belgium: Von Karman Institutdefor 

Fluid Dynamics, 1965) VKI Course Note 55b, pp. G46-50, 
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For Mach numbers of 0,8 or less, the discharge angle can be 
calculated quite accurately by 



T 3 -*- 

L 

where: 



a ■ minimum throat width 
t * blade thickness 

y 



then with r f - and using Eq. 140 in Eq. 139, 

- -a 

€ i J f )] 



which can be expressed as 



v t , \ ~ £-z 



1 



Assuming constant pressure at the throat, 

- - 3l 

f ' p t 



and 



where 



Ah ~ T 

_ -r 

Tth. . 'to _ Jto z °‘l c — 






In Eq. 144 the term X is 



Tt. 



r to ; 



*-l 



= i-3k = i-fJkV 



The boundary layer is assumed to be turbulent and the profile can 
therefore be expressed by 



’th 



•WV 



(140) 

(141) 

(142) 

(143) 

(144) 

(145) 

(146) 
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Then with the displacement thickness of the boundary layer defined as 






- (l-X, 



I — Xfi 



(147) 



Eq. 142 can be written 



M. 






(148) 



The loss coefficient to the throat is 



<♦*1 

j v * 



(149) 



th 

The loss coefficient can also be expressed in terms of the kinetic 
energy lost > by 






C4, 4 £ 

* Z rrv 

where E represents energy rate at the discharge due to the average 

velocity V^, or 

Ey«^[ a -^]T + £/(“/-r)ty 



(150) 



(151) 



Equation 151 can be rewritten to give 




Then using the energy thickness which is defined as 




(152) 



(153) 



(154) 
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The restriction factor J represents that part of the throat 
opening in which would occur the uniform theoretical velocity, 
therefore 




Now defining an energy parameter H*** by 



,***_ ^ 



Eq. 155 becomes 

H*** 



- I 



5 “ 

The loss coefficient Jjr in Eq. 157 accounts for the losses that 
occur from the inlet to the throat of the blade channel. No means 
exists for predicting this loss coefficient. Half the total loss 
coefficient for the blade row 



losses. Eq. 157 then becomes 



will be used to represent these 



l- 



H** - \ 

2L 



,13 



4. Method of Evaluating H*** 

By toC of the binomial theorem, the denominator of the integral 
part of Eqs , 147 and 153 can be expanded f yielding 

v-1 



+ A e 



y Ac y 



ft sj Zi ^\ x vi 2*n. v/ a 

(j-Xey j - \ y y 

The integral of Eq. 153 is then 

ft 3™-. ^ 7^ .3 ^ / Hk 1 I. 

J + Xey Xe y + Xe. y +-Xe y + J ^ f 



Then, integrating and evaluating gives 

'' i I Xe Xe £ 



( jl 2 j = _L_ 



Xe 3 



J'-M 

Equation 153 can now be expressed as 



Xe 

n^ + i 



+ Jk. + 1*L . 

* l 5^+1 7*^>l \lw\+l / 
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Eckert, op. cit . , pp. 



159-160. 



(155) 

(156) 

(157) 

(158) 
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or, 




In a similar manner, the expression for Eq. 147 is 



sljt -|V-L + J_ U hi I ** , X e 4 , ^ 

£ e /yx e - 1 **+ 1 + 3 m+1 Srvv+I 7*t + l / 

The energy parameter can now be expressed by 



J ^ J. + Xg + X e ^ Xa. ^ 

****_ Xa~l 3*\ + | 7«t+-l il»v\ + l 

J— * -L_ v J<*. + 

X e ”l ** +1 3m+l 7**t + l e^v+1 



(159) 



The number of terms in the numerator and denominator of Eq. 159 
is considered sufficient to give good convergence for H***. 

Although the exponent m for the turbulent boundary layer is 
dependent on Reynolds number, it is taken to be a constant for 
this analysis. 
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APPENDIX B 

COMPUTATION OF OUTLET ANGLES AND LOSS COEFFICIENTS 



1. Outlet Angles 

The absolute discharge angles for the stator and the relative 
discharge angles for the rotor are computed by using the same 
methods. For these calculations, the assumption is made that the 
outlet angles are not influenced by the flow incidence angles. 

The factors which do affect the discharge angles are: 

1. Blade geometry; this results in the angles being a func- 
tion of the radius since spacing and possibly profiles change from 
hub to tip. 

2. Radial tip clearance; the effects of tip clearance are 
assumed to influence the rotor flow from the mean streamline out 
to the tip > with the largest effect being near the tip. 

3. Exit Mach number; in accordance with the experimental 
results surveyed by Ainley, the outlet angles are a function of 
Mach number. Values of the angles are calculated for Mach numbers 
M of 0.5 and 1.0. Then a smooth curve of outlet angle versus M 

is drawn between these points with an inflection point at M*0.75. 
Below h-0.5, the flow angles are assumed to be equal to the value 
computed for M»0.5. 

Vavra’s formula is used for the first approximation of the 
stator outlet angles for M-0.5. 




( 160 ) 



where 




(161) 



The effects of blade curvature are next taken into account 
by using the method given by Ainley,^' 




(162) 



14 



Ainley, op. cit . , pp. 3-4 
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II f< 

In Eq. 162, e is the mean radiu,' of curvature of the upper surface 
of the blade profile between the passage throat and the trailing 
edge. This quantity is approximated by 




where j and b are shown in Fig. 46. 




The discharge angle for M«1.0 is found in a manner similar 
to that described by Ainley, 

cos -t*) 

The change in OC due to the increase in local Mach number from 
0.5 to 1.0 is assumed to be constant along the blade height and 
equal to that computed for the mean radius. 

5 0C T - S 0( H = ~ 

Relative discharge angles for the hub and mean radii of the 
rotor are found by the same method as used for the stator. The 
effects of tip clearance are accounted for in the computation of 
^ for the tip position. With Ainley’ s formula, the flow 
angle for t ^ ie tip * s 






(163) 



(164) 



(165) 



(166) 
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The new terms In Eq. 166 are* 
h - blade height 

discharge angle before tii clearance is accounted for 
X - factor which depends on the type of shroud 

Ainley suggests that X be set equal to 1.35 if no shroud is 
arranged at the rotor tips. However, Ainley' s value has been 
determined with a mean streamline analysis. To localize the 
effects of tip clearance, X has been changed to 5.4, a value which 
gives about the same overall tip clearance effect as if the value 
of X=l,35 were used for the mean streamline. Fig. 47 shows that 
the areas under the two lines representing X as a function of radius 
are approximately equal. The slope of the solid line representing 
X for this method is obtained because of the assumed linear varia- 
tion of outlet angle with radius between the mean radius and the tip. 




Tip Clearance Factor X as 
Function of Radius 

Fig. 47 

Sample calculations for stator outlet angles are shown in 
Table I. A sample of rotor outlet angle computations can be seen 
in Table II. 
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TABLE I 



SAMPLE CALCULATIONS 
MOD I turbine 


FOR STATOR OUTLET 
M x <: 0.5 


’ ANGLES 




radius - in. 


3.597 


4.125 


4.950 


spacing (s) - in. 


1.738 


1.994 


2.392 


throat opening (a) - in. 


0.455 


0.590 


0.804 


trailing edge thickness (t) 


- in. 


• 0.044 - 




a 

S 


0.2619 


0.2960 


0.3360 


0100 =A 


21.4 


13.6 


7.48 


K t -l-<0.0027) (A) (|) 


0.9849 


0.9891 


0.9932 


a 


sK t 


0.2658 


0.2991 


0.3384 


r,* -1. a . 

* ” cos ( iKT 5 


74.6° 


72.6° 


70.2 




0.703 


0.9106 


1.37 


a. = oc*+ 4(|) 


75.3° 

Mj®!. 0 


73.5° 


71.6 


Use 4.125 in. from centerline as representative radius 


. 


Projected trailing edge thickness (t*)-in. 
- cos -1 ** deg. 

= 73.5 - 71.5 = 2.0 


0.13 
71. 5* 




oC 


e 

73.3 


71*5° 


69.6° 
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ABLE II 



SAMPLE CALCULATIONS FOR 
MOD II turbine 



ROTOR OUTLET ANGLES 

tip clearance “0.033 in. 



M * 0.5 



radius - in. 




3.3235 


4.197 


4.918 


spacing (s) - in. 




1.595 


1.4643 


1.7158 


throat opening (a) - in. 
trailing edge thickness 


(t) - in. 


0.382 


0.560 

0.048 


0.721 






a 

s 




0.3295 


0.3824 


0.4202 






4.1397 


3.278 


2.7975 


(4— r=A 




109.0 


50.3 


29.8 


K t = 1- (0.°O2?)(fl)(-|-) 




0.9030 


0.9481 


0.9662 


a 




0.3649 


0.4033 


0.4349 












-68.5° 


-66.2° 


-64.2 






1,75 


2.37 


3.04 






. -70.3° 


-68.6° 


-67.2 


k - in. 








0.033 


k /h 








0.0207 


X*(4) (1.35) -type shroud 


factor 






5.4 


^ = 50° 








0.86603 


COS/^' ,/3 z '= -C7.Z° 








0.38268 










0.25297 


tan £ 








-2.4142 


(l-B)tan^' 








-1.8035 


tan ^ 








0.57735 


(1-BJtan/^B TflU^ = 








-1.6574 


/^ 2 «tan _1 (C) 








-58.9* 
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TABLE II (continued) 

SAMPLE CALCULATIONS FOR ROTOR ANGLE 

m 2 «i.o 



Use radius = 4.197 in. as representative radius. 
Projected trailing edge thickness (t*)-in. 

c<ps X^r) 

*-68.6+65.6= -3.0 



0.108 

-65.6° 



radius - in. 

A 



3.3235 4.197 4.918 

-67.3° -65.6° -62.6° 
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2. Loss Coefficients 

Loss coefficients were predicted using Ainley's methods. 

The techniques used are completely described in Refs. 1 and 2. 

The assumptions, definitions, and equations which are necessary 
for a basic understanding of the method are described below. 

It is assumed that loss coefficients are independent of Mach 
number and, for a given glade geometry, are a function of the flow 
incidence on the blade. Losses are divided into profile losses, 
secondary losses, and losses due to tip clearance. Mixing losses 
are not accounted for separately but are included in profile and 
secondary losses. 

A parameter Y is defined by 

Y = Loss in total pressure 

Total pressure at outlet - outlet static pressure 



For incompressible flow, which is assumed for these calculations, 
the loss coefficient is 



3 = 






I *■ £ Yji 

The subscript y on Y in Eq. 168 indicates the type loss repre- 
sented; i.e,, Yp would be the profile loss. 

For zero incidence, the Y for profile losses is 



Y P 












0 . 2 . 






where : 



^ - blade discharge angle not accounting for tip clearance 

t - maximum blade thickness 
m 

c - blade chord 

Equation 169 is applicable to both the stator and rotor. However, 
for stator use the blade inlet angle is zero. The quantities 

^ and Yp^.^are taken from Fig. 4 in Ref. 2 where these 
quantities are plotted as a function of blade solidity. The maximum 
value of t- m j c is set equal to 0.25 for thick blades. 

To determine profile losses for incidences other than zero. Fig 

Vp L 



6 in Ref. 2 shows 



versus 

iPttro} <S 



where i is defined as that in- 



cidence where the losses are twice the minimum losses. 



(167) 



(168) 



(169) 
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Secondary and tip clearance losses are computed by 



Y s +Y* 



h J L VcJ co,*^ 



(170) 



where: 

/ ^ n _- the mean of the gas inlet and exit angles of the blade row 
6 - factor representing type shroud, similar to X 

in Eq. 171: 

■ (annulus area at discharge) (cos ) 

* (annulus area at inlet) (cos ^ ) 

I.D. and O.D. - inner and outer diameters, respectively, at the 
blade exit 

Values of are obtained from Fig. 8 in Ref. 2. Ainley suggests 

a value of 0.5 for B for the type shroud used. With the same rea- 
soning as previously used for the effect of tip clearance on dis- 
charge angles, B is set equal to 2.0 for the losses at the tip and 
zero for the mean radius and hub locations. 

Sample calculations for the prediction of stator loss coeffi- 
cients are shown in Table III. Sample calculations for prediction 
of stalling incidence and rotor loss coefficients are shown in 
Table IV. 
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TaBLE III 



SAMPLE CALCULATIONS FOR STATOR LOSS COEFFICIENT 



MOD I turbine 
radius * 3.597 in. 
pitch (s) * 1.738 in. 

All figures referred to are in Ref. 2. 



inlet angle ( 0( o ) = 0 
outlet angle ( 0^ ) =* 75.3° 
chord (c) = 2.691 in. 









V/VcY**, 
o.z I 



oU- o 



= °' 647 



PitcyI _ 1738 
CHORD ** 2.631 



= 4S9 






^ f [(V.,)‘/o . 1)]. f[c— ■)•/(, , ^ 

A = ^Yo.0334) A = 0,0056 s) 

TAM''0-^) = 4.2.3° 

=2[tam <x. 6 -taM C iscos = 2 [- (taM - 7S.3 0 )] cos 6?. 5 



X . D . = £ » 6 4 ri, 

O , D, - 3.538 4/1. 



Vc 

Y< = 



= 3.S455 



[O ■ O OS6] [3. S45 8] *j> s ll s -3‘~| = o 



= 0.0451 



Y-, + Yp = o. 0 929 

CORRECTION FACTOR For TRHaiMfi ED6E THICK MESS , F=103 

(fL$ ■?) 

(j.05)(o.0929) = OO 949 



3 = r. 



Q'Q 349 
\-bO. 094-9 



= 0.0666 
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TABLE IV 



SAMPLE CALCULATIONS FOR STALLING INCIDENCE AND 
ROTOR LOSS COEFFICIENTS 



MOD II turbine 

radius ■ 4.918 

t ■ 0.728 in. 
max 



s * 1.7158 in. 
k * 0.033 in. 
c ■ 1.967 in. 



A - 30° 

A'- - 67 - 2 ° 

blade height(h) * 1.595 in. 



Figures referred to are in Ref. 2. /? 

W(V^)*[Vvo-V^)' 
t'tiM- < - f: ■°- 444 ' 

Yp^^O.0365 +(_O,l9S3)(0.lS4- 0.f56S^Q.2^)°' 44<£, _ o.O(o(>Z 

ft ‘ * 

= 0.955 (f *.T) , £a., 



(J' a7S ) ^(|=0.7S) 

‘•(i. MS ) sK ' 5 ‘ ■ ^s = -V 0 (Fig.7) , l'^3?.S°-7°= 2S_S° 

Y* * Y* = [?> + B (-£)] [y%] 5^-o.oisei 

COS^‘ = 0.387S2 cost's 0,15017 

x_/J7/3 / \ 2 /, x D sl T.d. outlET= 3,323S crt, OP. outlet - 4.3t6 tA, 

' [\ /^V/V + cTcT)] X.D, IWLET = 5,324S {A. OP.lMl.ET = 4,8SO€ tV. 

41= 70, 3* 

Ik _ W4,9l8 a - S.3235 2 ) Co*7°3 

fl, tr(4,8sos T --5.32« I ) c. 0& 3o° 

7) = ^0,148) r h-o.ooiz. (F ;<j, 8) 

Y, = [0.0072 +O ,O,08l][ C / S/ J Z [^iy 

Y* + Y< = (_o, oo 3606) [VvJ* 



- oAZC 
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TABLE IV (continued) 



SAMPLE CALCULATIONS 


FOR ROTOR LOSS 


COEFFICIENT 




ft* 


-1.5 


0.0 


1.5 




-38.25 


0 


38.25 


/&>- L + fio 


-8.25 


30 


68.25 




-0. 1450 


0.5774 


2.5065 


TflN/S, _ a 

2L M 


-1.2619 


-0.9008 


0.0638 


/C= rm' 1 (p) 


-51.6 


-42.0 


3.6 




0.6212 


0.7430 


0.9980 


TAM^-TArt/^'* B 


2.2339 


2.9563 


4.8854 


%c = 2 (B)c 0 s/^ 


2.7752 


4.3928 


9.7510 




7.7016 


19.2964 


95.0813 


Vcos’/U 2 


4.1726 


2.4385 


1.0061 




0.2345 


0.3434 


0.6980 


VY Pt=c > (**<•> 


2.1 


1.0 


4.5 


Y P 1 IV%-J Y Pt--») 


0.1390 


0.0662 


0.2979 


Vp’-Ys +Y k 


0.3735 


0.4096 


0.9959 


(CORRECTION F/)CTO?t)(£ V)=- C ' 


0.3922 


0.4301 


1.0457 


C 

' 1 + c 


0.2817 


0.3007 


0.5112 
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APPENDIX C 



COMPUTER PROGRAM 

A computer program using Fortran IV was written for use with 
the performance analysis „ The program was used to predict perfor- 
mance values for the MOD I and MOD II turbines. However, an 
attempt was made to keep the program general. If the methods of 
this thesis are utilized for expressing turbine characteristics, 
this program can be used for other single stage axial turbines, 

Input for the program consists of information representing 
the turbine characteristics, an indicator specifying the detail 
of output desired, and the conditions for which performance values 
are to be obtained. 

Twenty input cards are used to introduce the turbine 
characteristics. The information contained on each of these cards, 
with required dimensions, is listed below. 

1. Stator mean streamline gas outlet angles (ALFAM- radians) 
for mean streamline exit Mach numbers of 0.5, 0.7, 0.75, 0.8 and 
1 . 0 . 

2. Relative rotor mean streamline gas outlet angles (BETAM- 
radians) for the same Mach numbers listed in (1.). 

3. Radii ahead of the stator (RC-in.) for the five stream- 
lines such that the flow area is divided into four equal parts. 

4. Assumed radii of the five streamlines at the stator exit 
(RS-in.). 

5. Same as (4.), only for the rotor exit plane (RR-in.). 

6. The predicted stator loss coefficients (ZETAS) for the 
hub, mean radius, and tip. 

7. Ten values of incidence ratio (RINC) ranging from -2,0 
to 1.6 jin increments of 0.4. 

8. Ten values of rotor loss coefficients (ZETAR1) for the 
hub corresponding to the incidence ratios of (7.). 

9. Same as (8.), only for the mean radius (ZETAR3). 

10. Same as (8.), only for the tip (ZETAR5) . 



112 



11. Length L shown in Fig. 2 (CL- in.) and the curvature 
factor K (CK) . 

12. Differences in stator gas outlet angles (DALF-radians) 
from the angle for the mean streamline, for streamlines 1, 3, and 5. 

13. Same as (12.), only for the relative rotor outlet angles 
(DBET) . 

14. Ten radii for the stator exit plane, equally spaced or 
approximately equally spaced, ranging from the radius of the hub 
to the radius of the tip (Rl-in.). 

15. Ten values of stator throat opening (Al-in.) corresponding 
to the radii of (14.). 

16. Same as (14.), only for the rotor exit plane (R2-in.). 

17. Ten values of rotor throat opening (A2-in.) corresponding 
to the radii of (16.). 

18. Inlet blade angles for the rotor (BETO-degrees) for the 
hub, mean radius, and tip. 

19. Ten stall incidence angles (STALI -degrees) corresponding 
to the radii in (14.). 

20. The radial tip clearance of the rotor (TIPC-in.), the number 
of stator blades (ZS), and the number of rotor blades (ZR) . 

Input card number 21 specifies the output to be printed, and its 
use will be described later. The remaining input specifies the con- 
ditions for which performance values are to be found, and enters the 

estimations of the flow Mach numbers used for the first approxima- 
tions in the iteration process. Input card number 22 specifies the 
number of sets (NSETS) of operating conditions for which solutions 
are to be found. A card containing the following information is 
used for each point specified by NSETS: 

a. Estimated Mach number of the flow ahead of the stator (AMC). 

b. Estimated Mach number of the mean streamline flow after the 

stator (AMS). 

c. Total inlet pressure (PTO-psi.). 

d. Total inlet temperature (TT0-°R) . 
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e. Estimated relative Mach number of the mean streamline 
flow at the rotor exit (AMR) . 

f. The rotor speed (RFM) . 

g. The ratio of total inlet to static discharge pressure (PR). 

There are eight subroutines in addition to the main or execu- 
tive part of the program. The subroutines and their main functions 
are listed below: 

1. Subroutine PARAB is used to determine the coefficients of 
the parabolic equations used to approximate curves. 

2. Subroutine LSQPOL determines by the method of least squares 
the coefficients of the fourth order polynomials used to approxi- 
mate curves. 

3. Subroutine CHAN computes the flowrate and the reference 
flowrate from conditions ahead of the stator. 

4. Subroutine STATOR determines the axial velocity ratios that 
satisfy the equation of motion at station 1. 

5. Subroutine FLQWREF computes the reference flowrate at the 
stator exit and the rotor exit^and adjusts the axial velocity of the 
mean streamline flow at these locations to satisfy overall continuity, 

6. Subroutine SLINE checks streamline continuity for the stator 
exit and the rotor exit*and determines new streamline radii to satisfy 
streamline continuity. 

7. Subroutine R0T0R1 converts the absolute flow properties 
ahead of the rotor to relative flow properties. 

8. Subroutine R0T0R2 determines the axial velocity ratios at 
station 2 that satisfy the equation of motion. 

The main program computes the mass- flow-weighted value of total 
inlet to static discharge pressure ratio, and adjusts the Mach number 
of the flow ahead of the stator to obtain the pressure ratio speci- 
fied. The overall turbine performance values of ^ , M^,H.P.,r*, and 
k^ s are also computed in the main program. 

Computer output representing the solution for each set of condi- 
tions is printed on two pages. The flow properties at the stator 
exit plane are printed on one page . T A sample of this output can 
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be seen in Table V. The second page contains the flow properties 
for the rotor as well as the overall turbine performance values. 

A sample of the second page of output is in Table VI. Host of the 
output quantities are self-explanatory. However, the following 
symbols may not be readily recognized and are defined as: 

SLINE - streamline 

ZETAP - one-half the total loss coefficient, used for calcula- 
tmg | 

A - throat opening of the blade channel 

W-FRAC - fraction of the flowrate w between the hub and the 
indicated streamline 

BO - rotor blade inlet angle 

INCID - rotor flow incidence angle 

W 

PSIR - relative flow velocity coefficient 2 

W 2th 

In addition to the main output there are additional WRITE 
statements in the program. The additional output allows the user 
to follow more closely the intermediate steps of the solution pro- 
cess. This output is also helpful when program changes are made 
which require debugging. The extra output is not desirable for pro 
duct ion runs because of the large amounts of computer running time 
and printout that result. When the additional output is wanted, 1 
is read into the computer on input card number 21 for the indicator 
IND. If 0 is read into the computer for IND, only the main output 
representing the final solution will be printed. 

The statements of the program are listed on the following pages 
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